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We discuss the spectral, transport and magnetic properties of quantum nanowires composed of N ^ 13 
atoms and containing either even or odd numbers of valence electrons. In our approach we combine Exact 
Diagonalization and Ab Initio calculations (EDABI method). The analysis is performed as a function of 
the interatomic distance. The momentum distribution differs drastically for those obtained for even N with 
those for odd TV, whereas the Drude weight evolve smoothly. A role of boundary conditions is stressed. 
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For low-dimensional systems, the procedure starting from the single-particle picture (band structure) 
and including subsequently the interaction via a local potential may not be appropriate. In this situation, 
one resorts to parametrized models of correlated electrons, where the single-particle and the interaction- 
induced aspects of the electronic states are treated on equal footing. The single-particle wave-functions 
are contained in the formal expressions for model parameters and should be calculated separately. We have 
proposed JQ to combine the two efforts in an exact manner. 

A separate question concerns the role of boundary conditions (BCs) in atomic rings, particularly under 
the presence of spin frustration for odd number of atoms. This problem was investigated numerically Q 
and the optimal BCs for a correlated system were found to remain usually the same as for the ideal Fermi 
gas on the lattice. However, the general proof of this basic fact have been elaborated very recently 1 3 1. 

In our method of approach (EDABI), we determine first rigorously the energy of interacting particles 
in terms of the microscopic parameters for a given BCs and then allow the single-particle wave functions 
(contained in the parameters) to relax in the correlated state. The method has been overviewed in 1 4 1, so we 
concentrate here on its application to nanochains of N ^ 13 atoms, containing either even or odd number 
of electrons. The discussion of parity effects complements our recent study of correlated nanochains |5 1, 
where we consider the properties of even-iV systems only. 

We consider the system of N lattice sites, each containing a single valence orbital and (i.e. hydrogenic- 
like atoms). Including all long-range Coulomb interaction and neglecting other terms, one can write down 
the system Hamiltonian in the form 

H = ef n i + * ( e ^ HNc l^+^ + h - c ) + Uj2 n ^ n n + J2 K ij Sn i^j, (1) 

j ja i i<j 

where 5ni = rii — 1, e° ff — e a + N^ 1 J2i<j(2/Rij + Kij) ( m Ry) 1S the effective atomic level, Rij 
is the distance between the i-th and j-th atoms, t is the nearest-neighbor hopping, cf) is the fictitious 
(dimensionless) flux through the ring, U and Kij are the intra- and inter-site Coulomb repulsions. The 
last term represents the correlated part of the long-range interaction. One can easily show, that the unitary 
transformation cj a — > e~ 1 ^ 'l N ' Cj a allows to accumulate all the complex phase factors in the terminal 
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Fig. 1 Parity effect on spin ordering: spin-spin correlations for nanochains of N = 10 (a) and TV 
The values of the interatomic distance R are specified in the atomic units (ao = 0.529 A). 



11 (6) atoms. 



hopping term [6 1, which than takes the form t(e~ l ^c\ a CNa + h.c.) and can be regarded as generalized BC 
Such form is particularly convenient for numerical purposes, since majority of the hopping terms are real. 
Hereinafter, we do not distinguish between the system with a fictitious flux and with generalized BCs. 

The Hamiltonian Q is diagonalized in the Fock space with the help of Davidson technique Q, As the 
microscopic parameters e~, t, U, and are calculated in the Gaussian basis composing the Wannier 
functions, the orbital size of the ls-like state is subsequently adjusted to obtain the minimal ground state 
energy Eq as a function of the interatomic distance R. 

We now discuss the spin-spin correlations (Si • Sj) in the system ground state as a function of the 
discrete neighbor distance and R, as presented in Fig.^ The effect of spin frustration in the half- 

filled case (TV e = TV) is remarkable for large R, where the quasi-long range order for TV = 10 (cf. Fig. 
Qt), indicating the power-law decay of (Sj • Sj) for the Heisenberg spin chain, disappears for TV = 11 (cf. 
Fig.[IJ)), where we observe a fast, exponential decay. For small values of R, however, the effect is weaker, 
since the antiferomagnetic order is reduced by charge fluctuations [5 1. We also observe, that the values of 
the spin gap (not shown) are significantly higher for odd TV in the large- R range, what can be explained by 
the fact, that the ground-state energy of the Heisenberg antiferromagnet is of the order Eq ~ JS(S +1), 
where J = At 2 / (U — K) is the kinetic-exchange coupling parameter and S is the total spin value. One 

can expect now, that the following inequality is satisfied E G 2 — E G ~ 2 > E G =1 — E G =0 , where the left- 
and the right-hand sides represent the spin gap for the odd— and the even—N systems, respectively (both at 
the minimal-spin configuration). The detailed behavior of the system spin, as well as the charge and the 
optical gaps will be discussed elsewhere. In the remaining part of this paper we focus on the parity effect 
for the Fermi-Dirac distribution function and the transport properties. 

The electron momentum distribution for nanochains with optimal BCs is shown in Fig.|2k. The discrete 
momenta, corresponding to the solutions of the single-particle part of the Hamiltonian Q for a finite N, 
are given by 

M^) = ^T^ 0^q<N. (2) 

The optimal BCs, corresponding to the minimal ground-state energy Eq, are realized for = when 
N = An + 2 (periodic BC), <fi = tt when N — An (antiperiodic BC) and <fi — tt/2, 3tt/2 when is odd 
l3l . A basic analysis of Eq. (0 shows, that for the optimal BCs, the Fermi momentum state kp = ir/2R 
is newer reached for even TV, whereas for odd TV it happens, for a single value of q. This circumstance 
has tremendous implications for electronic structure of a nanochain, however, of almost does not effect its 
transport properties, as we show in the end part of this paper. 
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Fig. 2 Electron momentum distribution for chains of Fig. 3 Normalized Drude weight for nanochains in 
N = 8 -5- 13 atoms: (a) datapoints for optimal boundary the half-filled case (a), and for a system with a sin- 
conditions (BCs) and the sample n(k) curves averaged gle hole (6). The optimal boundary conditions are 
over BCs (solid lines); (b) the original n(k q (<j})) func- specified for each curve, 
tions used for the averaging. 



The remarkable feature of these results, is that the datapoints for different but even N (cf. full symbols 
in Fig. [2^) locate smoothly on the universal curve (for each R), when optimal BCs are applied. This is not 
the case for odd N (cf. open symbols in Fig.|2^), when the systematic dependence on N suggests quite slow 
convergence to the even-TV results with the increasing N. However, the discussion of N — > oo limit is 
beyond the scope of this paper, since we concentrate here on nanochains. To analyze such systems in detail, 
we have displayed in Fig. |2j) the continuous electron momentum distribution obtained for the dense set of 
k g ((f)) defined by Eq. (|2ji, when <\> G (0, 2tt). The datapoints corresponding to optimal BC for N = 10 and 
11 are also presented to show, they are situated close to the different local extrema of n(k q (4>)) (e.g. maxima 
for even N and minima for odd, or vice versa). Except of different frequency of internal oscillations (equal 
to 2NR), the n(k q (<f))) functions for N — 10 and 11 looks almost identically. This is the reason, why 
various physical properties of small clusters are often averaged over BC, particularly in 2D [6|. We also 
perform such averaging to obtain almost size-independent n(k) functions, drawn again in Fig.|2k (the 
details of the averaging procedure will be published elsewhere). However, the elimination of the internal 
oscillations for a given N may only be considered as an approximation of N — > oo scaling procedure, and 
in the case of momentum distribution n(k) seems less accurate then fitting the Luttinger-liquid formulas to 
even N data, which we proposed before [5 1. Nevertheless, the common nature of either original n(k q {<p)) 
or averaged n(k) functions for both even and odd N, illustrated in Fig. |2] helps to understand why the 
chain parity does not effect its Drude weight even for optimal boundary conditions, when the structure of 
the momentum space is significantly different. 
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The normalized Drude weight, shown in Fig. [5] is defined in the standard manner | 

1 d 2 E G 



D 



(T) dcj) 2 



(3) 



where (T) is the average kinetic energy and <^> m in denotes optimal BGs. In the half-filled case N e = N (cf. 
Fig-Et) Drude weight gradually decrease with N, as we have shown for even N Q. The most interesting 
feature of these results is, that the curves for odd N fits smoothly between those for even N, with very weak 
parity effect (totally incomparable with that present in the charge gap and electron momentum distribution, 
when optimal BC are applied). This observation can be understood when we take into account, that the 
Drude weight defined by Eq. |3]is the integral quantity, involving the summation over all the excited states 
of the Hamiltonian Eq. Q, so it cannot be determined only by the electronic structure near the Fermi 
points, particularly for a small system. 

The parity effect on Drude weight disappears for the system with a single hole (N e = N — 1, cf. Fig. 
fJJ)), in which the magnetic frustration is absent. For this case, the Drude weight evolution with R is very 
interesting. In the weak-correlation range (R/ao < 2) the chain shows a highly-conducting behavior 
for each N. Next, in the intermediate range {R/ao = 4 -j- 5) the Drude weight decrease rapidly with 
N, indicating an insulating (Mott-Hubbard) state in the large- N limit. In the strongly-correlated range 
(R/ao ~ 10) the Drude weight approaches again its maximal value D = 1. Such a behavior can be 
explained, when we analyze the situation in two steps: First, for low values of R, the bandwidth-to- 
interaction ratio is small, and the system with a single hole does not differ significantly from a half-filled 
one. This is why in both cases Drude weight decreases gradually with both N and R, as the tunneling 
amplitude through the barrier of a finite width. Second, for the largest values of R, the system can be 
described by an effective t — J model |9| with a coupling constant J = 4t 2 /(U — K) <C \t\ (where 
K = Kj t j+i denotes the nearest neighbor Coulomb repulsion), which corresponds to an asymptotically- 
free hole motion. Then, it become clear that in the intermediate range the Drude weight has to suppressed, 
what can be interpreted in terms of a partially localized spin-ordered state. It would be very interesting to 
test experimentally this result, possibly for a mesoscopic atomic ring. 

In summary, we have shown that a nanochain parity effect strongly its electronic structure and mo- 
mentum distribution, and that the effects are opposite for these two principal characteristics. Namely, the 
presence of a discrete momenta in the Fermi point reduce significantly the finite-size effects on the system 
charge-gap, but amplifies them in the case of momentum distribution. On the other hand, the parity effect 
is weak, or even absent in the case if system transport properties. Additionally, an interesting crossover 
behavior has been identified for the chain with a single hole, for which the quantum-liquid regions are 
separated by a partly localized state. 
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